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Abstract 

Isotropization occurs on time scales much shorter than the thermal 
equilibration time. This is a crucial ingredient for the understanding 
of collision experiments of heavy nuclei or other nonequilibrium 
phenomena in complex many body systems. We discuss in detail 
the limitations of estimates based on standard "linear" or relaxation- 
time approximations, where isotropization and thermal equilibration 
rates agree. For a weak-coupling <^ 4 -model the relaxation-time 
approximation underestimates the thermal equilibration time by 
orders of magnitude, in contrast to the isotropization time. The 
characteristic nonequilibrium isotropization rate can be enhanced as 
compared to the close-to-equilibrium value. Our results are obtained 
from the two-particle irreducible effective action, which includes off- 
shell and memory effects and does not involve a gradient expansion. 
This allows us to determine the range of validity of a description 
to lowest-order in gradients, which is typically employed in kinetic 
equations. 
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1 Introduction and overview 



Understanding the dynamics of quantum fields far away from the ground 
state or thermal equilibrium is a challenge touching many aspects of physics, 
ranging from early cosmology or collision experiments with heavy nuclei 
to ultracold quantum gases in the laboratory One of the most crucial 
aspects concerns the characteristic time scales on which thermal equilibrium 
is approached. Much of the recent interest derives from observations in 
collision experiments of heavy nuclei at RHIC. The experiments seem to 
indicate the early validity of hydrodynamics, whereas the present theoretical 
understanding of QCD suggests a longer thermal equilibration time. 

However, different quantities effectively thermalize on different time 
scales and a complete thermalization of all quantities may not be necessary 
to explain the observations. We have pointed out in Ref. £Q that 
the prethermalization of crucial observables may occur on time scales 
dramatically shorter than the thermal equilibration time. In particular, an 
approximately thermal equation of state may be reached after an extremely 
short time. From this early time on the effective kinetic temperature 
is already very close to the final equilibrium value of the temperature 
and the equilibrium relations between average pressure, energy density 
and temperature hold. Beyond the (average) equation of state, a crucial 
ingredient for the applicability of hydrodynamics for collision experiments j2] 
is the approximate isotropy of the local pressure. More precisely, the diagonal 
(space-like) components of the local energy-momentum tensor have to be 
approximately equal. Of particular importance is the possible isotropization 
far from equilibrium. The relevant time scale for the early validity of 
hydrodynamics could then be set by the isotropization time 1 El • 

We argue in this paper that isotropization generically happens on time 
scales much shorter than the thermal equilibration time. 2 For this we 

lr The use of hydrodynamics is based on conservation laws, which involve the different 
components of the energy-momentum tensor. If dissipative terms can be neglected, these 
equations can be closed if an isotropic pressure can be expressed as a function of the energy 
density by an equation of state. We have shown in a previous publication 1 that the ratio 
between pressure and energy density settles already at a very early prethermalization time. 
Therefore, isotropization sets the relevant time scale for the use of hydrodynamics, in the 
absence of large dissipative effects as suggested by current experiments at RHIC 0. 

2 A possible mechanism for fast isotropization in QCD in terms of plasma instabilities 
in high-temperature gauge theories has been proposed in Ref. [3J- We do not consider 
the interesting topic of plasma instabilities here. For discussions of instabilities in a QCD 
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concentrate on a general framework for the description of isotropization 
in far-from-equilibrium quantum field theory jjj. Already for a simple 
scalar theory with small quartic self-interaction we find that there is a 
large separation of times between isotropization and thermalization. We 
compare our results with the most common estimates, which evaluate a 
relaxation time for small excitations away from equilibrium. This so-called 
"linear" or relaxation-time approximation jS] assumes that there is a more 
or less universal characteristic rate for many out-of-equilibrium properties. 
In particular, in the relaxation-time approximation the characteristic rates 
for isotropization and thermalization agree. In contrast, we find that for 
rather general initial conditions the relaxation-time ansatz fails even to give a 
correct order-of-magnitude estimate of the thermalization time r eq on which 
the Bose-Einstein distribution is reached. However, it describes well the 
(on-shell) isotropization rate, 7j~ , close to equilibrium at sufficiently late 
times. For earlier times, we find that the characteristic nonequilibrium 
isotropization rate can exceed the close-to-equilibrium rate for large initial 
anisotropies. We go beyond the relaxation-time approximation by calculating 
the nonequilibrium dynamics from the two-particle irreducible (2PI) effective 
action to three-loop order which includes scattering, off-shell and memory 
effects. The nonequilibrium evolution is solved numerically without further 
assumptions. In particular, we do not apply a gradient expansion. This 
allows us to determine the range of validity of a description to lowest-order in 
gradients, which is typically employed in kinetic equations and to motivate 
the relaxation-time ansatz close to equilibrium. The earliest time for the 
applicability of the lowest-order gradient expansion is found to be of the 
order of the characteristic isotropization time r; so . Many of the general 
findings such as a large separation of time scales between isotropization 
and thermalization, as well as the shortcomings of the relaxation-time 
approximation to describe the linearized dynamics will also be present 
in more complex theories like QED or QCD. However, in gauge theories 
particular features such as plasma instabilities will strongly affect the relevant 
time scale for isotropization [3]. 

We consider a class of anisotropic initial conditions with a high occupation 
number of modes moving in a narrow momentum range around the "beam 
direction" p 3 = ±p ts , reminiscent of some aspects of colliding wave packets 
in the central collision region (cf. Sec. Ej). We compute the nonequilibrium 

plasma see Refs. Other approaches to address fast thermalization include Refs. 0. 
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Figure 1: Isotropization and thermalization. The lower curve shows the 
characteristic nonequilibrium rate for isotropization as a function of time t in 
units of the renormalized thermal mass tjir. At late times it is well approximated 
by the thermal equilibrium on-shell ratio t^ 1 = — $]g 0q ^/2u;( eq - , indicated by 
the square symbol. This verifies standard relaxation-time approximations for 
isotropization at sufficiently late times. In contrast to isotropization, the approach 
to thermal equilibrium is not well described by the relaxation-time approximation. 
As an example, the upper curve shows — T, g (t;u,pt B )/2uj, which characterizes the 
nonequilibrium evolution of the imaginary-part of the on-shell self-energy. Only 
at very late times, t ^> Ti so — W0/itir, it comes rather close to its thermal value. 
The dotted-line represents a late-time fit explained in Sec. El 

isotropization rate for a rather weak coupling of the g 2 4 -interaction. (For 
our figures we use g = 1/2 and units in terms of the renormalized thermal 
mass rriR.) In Fig. ^ we compare this rate with the relaxation-time 
approximation and with the thermal equilibration time. The lower curve of 
the figure shows the characteristic nonequilibrium isotropization rate for the 
momentum p ts as a function of time. For times t > 500/m^ the isotropization 
rate is well described by the standard relaxation-time approximation, i.e. by 
Tj^ 1 = —Y^^ / 2uj( cq \ where the imaginary part of the self-energy — Eg eq ^/2 is 
evaluated for on-shell u/ eq ) for momentum p ts . The thermal equilibrium value 
— Sg 6q ' ) /2u;( eq ) is indicated by the full square in Fig. ^ At early times, one 
observes that the relaxation-time rate can underestimate the characteristic 
nonequilibrium isotropization rate. 
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For comparison, the upper curve of Fig. ^ shows the nonequilibrium 
evolution of the imaginary-part of the self-energy or, more precisely, the 
on-shell ratio — E e (t; u},p ts )/2u (cf. Sec. EJ). Only at very late times this 
rate comes rather close to its thermal equilibrium value — £[> ec ^/2u/ eq ). This 
reflects the slow approach of the momentum distribution to the Bose- 
Einstein distribution, which will be discussed in Sec. One observes that 
thermalization is an extremely slow process, and to verify a closer approach 
in an actual numerical calculation would require to go far beyond the times 
shown in the figure. Our results clearly demonstrate that the relaxation-time 
rate — /2u^ cannot be used at all to estimate the correct thermalization 
time. Most importantly, they show that isotropization happens on time scales 
much shorter than the thermalization time. This is analysed in detail in this 
paper. 

The paper is organized a follows. In Sec. El we derive the nonequilibrium 
evolution equations from the 2PI effective action and discuss the class 
of anisotropic initial conditions used in the present work. We solve the 
approximation of a three-loop 2PI effective action numerically, without 
further approximations. For further analytical discussion, we consider in 
Sees. H3 and 0] the linearized dynamics for small deviations from thermal 
equilibrium. For this we perform a gradient expansion to lowest order and 
establish its range of applicability by comparing to the "full" result including 
all orders of derivatives in Sec. E3 This is employed to derive the standard 
relexation-time approximation in Sec. 0] The discussion of the numerical 
results is presented in Sec. Eland we end with conclusions in Sec. El 

2 Nonequilibrium evolution equations 

Consider the single particle excitation which is obtained by adding a 
particle with momentum p to a system initially in equilibrium. For a 
weakly interacting system, which has long-lived single-particle excitations, 
this situation may be described in terms of a Boltzmann equation for the 
distribution function N p (t) for quasi-particles with momentum p and energy 
u) p (t). 3 For a small perturbation around thermal equilibrium SN p (t) = 

3 Without loss of generality for the argument we employ here spatially homogeneous 
ensembles. 
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N p (t) — Np eq ^ the relaxation-time approximation yields [1113 

6N p (t)=6N p (0)e-i® t , ( 2 .1) 

with the damping rate ^(p) determined by the imaginary part of the thermal 
equilibrium self-energy, and equilibrium on-shell frequency LV p ec ^ 

according to 

V(eq)/ (eq) 

7(p ^-^i!^- (2 - 2) 

This standard relaxation-time ansatz is widely applied to approximately 
describe the linearized dynamics around thermal equilibrium. 

The question about the range of applicability of the relaxation-time 
approximation has to be addressed in nonequilibrium quantum field theory. 
There are various assumptions leading to the above interpretation of the 
thermal equilibrium quantity (|2.2j) . In order to verify them we have to go 
beyond the standard gradient expansion employed in kinetic descriptions. In 
addition, no quasi-particle assumption should be employed a priori. We will 
therefore have to solve the nonequilibrium time evolution including off-shell 
and memory effects. 



2.1 2PI effective action and initial conditions 

Systematic approximations to describe far-from-equilibrium dynamics as well 
as thermalization fOJ El El Ell El El El E Ej from first principles can 
be efficiently based on the two-particle irreducible (2PI) effective action [TH] : 

T[ ( P,G} = SM + ^^G- 1 + 1 -TtG \<P)G+ T 2 [<p,G]. (2.3) 

To r 2 [0,G] only two-particle irreducible diagrams contribute, i.e. diagrams 
which do not become disconnected by opening two lines. To be specific, we 
consider a scalar quantum field theory with classical action 

SM = J (^<j>(x)d,cf>(x) - ^V(x) - g 2 <P\x)^ (2.4) 

in the symmetric phase with vanishing field expectation value, i.e. = 0. 
The classical inverse propagator is iGq (x,y) = 5 2 S c \[(p}/5(j)(x)5(j)(y)\ ( f >=0 . 
Summation over repeated indices is implied and we use the shorthand 
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notation J" = f c dx° f d 3 x with x = (x°, x) and C denoting a closed time 
path along the real axis starting at the initial time t = 0. 

The equation of motion for G is given by the stationarity condition 
ST/SG = 0, which from Eq. ()2.3|) reads 



G~\x, y) = G \x, y) - £(x, y; G) . (2.5) 
The self-energy E is related to T 2 [G] = T 2 [0 = 0, G] as 

In order to make the real and imaginary parts explicit we use the 
decomposition identity [TT) ITU] 

G(x, y) = F(x, y) - ~p(x, y) sgn c (x° - y°) , (2.7) 

where F(x, y) denotes the statistical two-point function and p(x, y) the 
spectral function. The spectral function encodes the equal-time commutation 
relations: 4 

P(x, V) U°=2/° = , d x op(x,y)\ x o =y0 = 5(x - y) . (2.8) 

To obtain a similar decomposition for the self-energy, we separate £ in a 
"local" and "nonlocal" part according to 

£(z, y- G) = -i£W (x; G)5(x - y) + S^(x, y; G) . (2.9) 

Since S^ ^ just corresponds to a space-time dependent mass-shift it is 
convenient for the following to introduce the notation 

M 2 (x;G) =m 2 + S (0) (x;G). (2.10) 

The imaginary part of the self-energy, — S p /2, is determined by 

^ nl \x,y) = Y. F (x,y) - % -?. p (x,y) sgn c (x° - y°) . (2.11) 



4 It is also directly related to the retarded propagator p(x, y)Q(x° — y°), or the advanced 
one — p(x, y)Q(y — x°). 
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Nonequilibrium dynamics requires the specification of an initial state. 
While the corresponding initial conditions for the spectral function are 
governed by the commutation relations ()2.8|h the statistical function F(x, y) 
and first derivatives at x° = y° = have to be specified. Here we consider 
systems described by Gaussian initial density matrices. This represents no 
approximation for the (non- Gaussian) dynamics for times t > 0, but just 
constrains the class of initial conditions. We consider a situation with an 
initially high occupation number of modes moving in a narrow momentum 
range around the "beam direction" ^3 = pii = ±p ts . The occupation numbers 
for modes with momenta perpendicular to this direction, p\ + p\ = p\, 
are small or vanishing. The situation is reminiscent of some aspects of the 
anisotropic initial stage in the central region of collision experiments of heavy 
nuclei. 5 More explicitely, we employ a class of initial conditions parametrized 
as 

T7(+ + -m n (p) + 1/2 

F{t 1 ,t 2 ;p)\ tl =t 2 =o = , (2.12) 

Up 

with d tl F(ti,0;p)\ tl=0 = 0, d tl d t2 F(t 1 ,t 2 ;p)\t 1 =t2=o = [™o(p) + V 2 ] w p for 
uo p = ^p\ + Pjj + Ml and M 2 = M 2 (t x = 0). The initial distribution 

function n (p) is peaked around the "tsunami" momentum p ts with amplitude 
A and width a: 

n (p) = A exp j-i [pi + (| p „| - PtB f] \ . (2.13) 

With these initial conditions, for our scalar model the most general statistical 
two-point function can be written as 

F(ti,t 2 ;p) = F(ti,t 2 ;p±,p\\) , (2.14) 

and equivalently for the spectral function, as well as for the self-energies 
£ir(ti,t 2 ;p) = Z F (ti,t 2 ;p±,p\\) and T, p (t 1 ,t 2 ;p) = Y; p (ti,t 2 ;p±,p\\). 

All momentum integrals appearing in this paper will be regularized on 
a lattice, with the momentum cutoff A chosen such that A ^> p ts . The 
renormalization of 2PI effective actions has been discussed in detail in 
Refs. j2U EE 1221 12H] • For the weak couplings employed here we follow the 
lines of Ref. ^7] to ensure that the relevant length scales are larger than 

5 Other interesting scenarios include "color-glass" -type initial conditions with 
distributions ~ exp(— \/pj~/Q s ) peaked around p% = with "saturation" momentum Q s . 
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the lattice spacing. We note that quantities such as the renormalized mass 
or renormalized damping rates can be directly inferred from the oscillation 
frequency and damping of F(ti,t 2 ;p = 0) or p(t\,t 2 ;p = 0) (see Eq. (jHHj) 
below) [7j. 

2.2 Time-evolution equations 

The exact time-evolution equations for known self-energies and S p are 
obtained from the stationarity condition of the effective action, or Eq. ()2.5|) 
by convolution with G. Using the decomposition indentities (|2.7j) and 
(12. lip , these are coupled differential equations for the spectral and statistical 
functions [HI DEI: 

[d 2 +p 2 + M 2 (t 1 )]F(t 1 ,t 2 ;$ = ~ r&t'^it^t'-^Fit'M-p) 

Jo 

+ / dt'£ F 0M';pM^2;p), 

Jo 

[dl+f + M 2 {t x )]p{t 1 ,t^P) = - / V t^p(t',t 2 ;p). (2.15) 

J to 



They are causal equations with integrals over the time history, starting from 
the time to = at which the initial conditions of Sec. 12.11 are specified. 
A characteristic anisotropy measure can be chosen as 

AF{t u t 2 ;q) = F{t 1 ,t 2 ;p ± = 0, Pll = q) - F{t 1 ,t 2] p ± = q, Pll = 0) , (2.16) 

which vanishes for the case of an isotropic correlator. Its exact evolution 
equation for known self-energies reads according to ()2.15|) : 

[dl+f + M 2 {t x )] AF(t 1 ,h;q) = 

- / dt' [S p (t 1 ,t / ;g)AF(t',t 2 ;g) + AS p (t 1 ,t';g)F(t / ,t 2 ;g)] 
Jo 

+ / dt' [E F (t x , t'- q)Ap{t', t 2 - q) + AE F (t u tf; qjp{t\ t 2 ; q)] . (2.17) 



Here we have defined the average 

F(ti,t 2 ;q) = ~ [F(ti,t 2 ;px = 0,p|| = q) + F{t u t 2l p± = q,p\\ = 



(2.18) 
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Figure 2: Diagrammatic representation of the contributions to ^[G] to three-loop 
order. 



and equivalently for the corresponding average of the spectral function, 
p(tx,t 2 ;q), and the corresponding average self-energies £ p (ti,t 2 ;g) and 

T lF (t 1 ,t 2 ;q). 

We solve Eq. ()2.15|) numerically using the coupling- or loop-expansion 
of the 2PI effective action to order g A or three-loop order without further 
approximations. To this order the effective mass term is 



M 2 (t 1 )=m 2 + 12g 2 f Ffat^P) 
Jp 

and the self-energies are JT] 



(2.19) 



F(tx,t 2 ;p- q-k) F{ti,t 2 ; q)F(t x ,t 2 - k) 



q.k 



-96g 
3 

-^p(ti,h;q)p(ti,t 2 ;k) 

-288s 4 / pih^p-q-k^Fih^qjFih^k) 



q.k 



1 
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p{ti,t 2 ,q)p{tiM;k) 



(2.20) 



using the notation f^=J d 3 p/(27r) 3 . The diagrammatic representation of 
the approximation is given in Figs. El and H3 We solve Eqs. (|2.15|) - (j2.20|) for 
the initial conditions numerically without further approximations along the 
lines of Refs. PHI [7j. 

In terms of these quantities the energy density e = T 00 and pressure 
components P l = T n [i = 1, 2, 3) are given by the diagonal elements of the 
energy-momentum tensor T^ v with 



dtidt 2 F{ti, t 2 ',q)\t!-- 



+ (q 2 + m 2 + 6g 2 F(t u h; p) F(t u h; q) 



[2.21) 





Figure 3: Diagrammatic representation of the corresponding contributions to the 
self-energy S. 



dt ' PAhA q)p(t', h] q) ~ Z p (ti,t'; q)F(t', t x] q)} , 
d tl d t2 F(tx, h- q) \ ti=h + 2q*F(t 1 , h; q) 

(2.22) 



- ( q 2 + m 2 + 6g 2 \ F(t u * l5 p) ) F(t u t x ; q) 



p 



+ 



~ jf jf 1 dt' [£ F (t 1; t'; #p(t', * l5 q) - Z p (t 1} t'; q)F(t>, t x] q)} . 



For an efficient (numerical) evaluation of these expressions it is advantageous 
to replace time integrals by time derivatives of F(ti,i 2 ;p) using Eq. ()2.15|) . 

2.3 Thermal equilibrium 

We finally note that out of equilibrium F and p are a priori not related 
by a fluctuation-dissipation relation. In contrast, in thermal equilibrium 
the periodicity ("KMS") condition for the propagator, G(x, y)\ x o =0 = 
G(x,y)\x°=-ij3, leads in Fourier space to the relation 6 

FM>(u,i?) = (n BE (co) + g^\u,p) , (2.23) 

with the Bose-Einstein distribution ^be(^) = l/[exp(/3o;) — 1] for inverse 
temperature f3. The equivalent relation holds for the statistical (Sf) and 
spectral (S^) part of the self-energy [TO]- We emphasize that we solve the 
nonequilibrium dynamics without any assumption about the validity of a 
fluctuation-dissipation relation. In order to compare the nonequilibrium late- 
time results to thermal equilibrium, we calculate in addition the thermal 
solution. This can be conveniently obtained from the p-part of Eq. ()2.15|) 

6 In order to have a real Q^ ec ^ (u> , p) we employ the Fourier-transformation as in Eq. 13.111 
below, which includes an additional factor of —i. 
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along the lines of Refs. (TSl C2] • The equation for the translation invariant 
thermal p^(s,p) with s = t — t' can be closed using Eq. (|2.23|) to eliminate 
with 

J dsF (eq) (s;p)cos(cus) = L E (w) + ^ dsp (eq) (s;p) sin(ws) . (2.24) 



3 Comparison with lowest-order gradient expansion 

For further analytical discussion — not to solve the equations — we consider 
the exact equations (j2.15|) in a gradient expansion 123 I2S EH E01 ED 
EH EH1 EH EH Hj to lowest order and use it to discuss the dynamics in 
the vicinity of thermal equilibrium. 7 This will also allow us to recover 
the relaxation-time approximation employing the corresponding additional 
assumptions described below in Sec. EJ Since the range of applicability of 
a gradient expansion is restricted to sufficiently homogeneous correlations 
at sufficiently late times, one cannot consider the early-time behavior in 
this case. We emphasize, however, that we are solving the nonequilibrium 
evolution using Eqs. (|2.15|) . which do not suffer from this restriction. In 
particular, our results can then be used to establish whether, and at what 
time, a gradient expansion can be applied. 

For sufficiently late times £i,t2 ^ the quantitative results turn out 
not to be notably affected by sending t — > — oo, which is required in a 
gradient expansion for practical purposes. Introducing relative and central 
coordinates, s = t 1 —t 2 and t = (ti+t 2 )/2, the Wigner transformed correlators 
are 

nt-,.,0, = /d.^F (t + ./2, t -./2;fl, 
g(t;co;p) = -i J dse iuJ3 p(t + s/2,t - s/2;p) , (3.1) 

where the factor i is introduced to obtain a real Wigner-space g(t; uj]p). Note 
that the time integral over s is bounded by ±2t because of the initial-value 
problem with ti,t 2 > 0. The equivalent transformation is done for the self- 
energies to obtain S^(t; o>, p) and T, g (t;u,p). In lowest-order of the expansion 
in derivatives of t and u, i.e. neglecting second derivatives, one finds from 

7 The following general analysis is not restricted to a coupling- or loop-expansion. 
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Figure 4: Range of validity of the gradient expansion. The solid line shows the 
on-shell S^AF + A£ e F - ^pAg - AHpg (thick line) and 2w^AF (thin line) in 
Wigner coordinates, as computed from the solution of the "full" equation (|2,17|) 
without a gradient expansion. If the gradient expansion to lowest order is correct, 
then both lines have to agree according to the AF- and A^-equations The 
inset shows the same results on a logarithmic scale to make the small 2uj^Aq 
visible. From these quantities we observe that the lowest-order gradient expansion 
becomes valid on a time scale similar to the isotropization time Ti so . 

the equations (j2~T5|) [2H ESI EE] : 

d 

2u—F(t;u,p) = E e (t;u,p)F(t;u,p) -i: F (t;u,p)g(t;u,p), 

2u^gfau,#) = 0. (3.2) 

We note that the structure of the gradient expanded equations is very 
similar to the exact equations (|2.15|) . In lowest order of the expansion the 
evolution equation for the spectral function g(t;u,p) becomes trivial. A 
frequently employed variant of the equations ()H.2j) also takes into account 
the change of the effective mass M 2 of Eq. ()2.19|) with respect to t. This 
contribution ~ (dM 2 / dt)(dF / du) and the Poisson brackets entering at 
second order, as well as all higher orders are neglected in the lowest-order 
expression (Q HI EE]. 

Following the same lines, in lowest-order of the gradient expansion one 
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finds from Eq. (|2.17|) the anisotropy equation 



E e (t; u, q) AF(t; u, q) + AE e (t; u, q) F(t] u, q) 
E F (t; w, q) Ag(t; u, q) - AE F (t; uj, q) g(t; u, q) , 
0. (3.3) 

In order to establish the range of applicability of the gradient expansion, we 
compare with the results obtained from Eqs. (j2.15j) - ()2.20j) which include all 
orders of gradients: We compute the Wigner-transforms of the results for 
AF(ti,t2',p) and Ap(ti,t 2 ;p) according to (j3.1|) . From this we then evaluate 
2lu^AF, and T> g AF + AT, g F — YlpAg — AHpQ, for on-shell uj. The results 
are represented as separate curves in Fig. 0] If the gradient expansion to 
lowest order is correct, then both lines have to agree according to Eqs. ()3.3|) . 
where these contributions appear as the LHS and the RHS of the equations. 
At early times the agreement is poor and the gradient expansion to lowest 
order cannot be applied as expected. We recall that because of the initial- 
value problem with t\,t2 > the time integral in the Wigner transform 
(13. lj) is finite. However, the effects of the necessarily finite time interval for 
a nonequilibrium evolution quickly become irrelevant and do not notably 
affect the curves. For times t > 100 /m R the two curves approach each 
other closely and the full result is indeed well approximated by the lowest- 
order contributions in the derivative expansion. This establishes the use of a 
gradient expansion for the anisotropy dynamics at sufficiently late times. We 
note that, in contrast to the differences AF and A^ measuring the anisotropy, 
the agreement was found to be much worse for F and g. Their description in 
terms of a lowest-order gradient expansion is very poor for the times displayed 
here. 

4 Linear and relaxation-time approximation 

Linearized dynamics around thermal equilibrium can provide a valid 
description at sufficiently late times, when the system is approaching thermal 
equilibrium. We denote the four- momentum p — (u,p) and write 

F(t;p) =F^\p) +5F(t;p). (4.1) 



2LU—AF(t;u,q) 



2uj^Ag(t;u,q) 
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Equation (|3.2j) can be linearized in the deviations SF(t;p) from the 
t-independent thermal equilibrium value F^ eq \p): 



d_ 

at 



5F(t;p) = - f S(p,q)5F(t;q) 

J Q 



r y, (eq) (r>) 

= - J S'{p, q) 5F(t- q) + 5F(t- p) , (4.2) 

with f q = J d 4 g/(27r) 4 . Here the "stability matrix" S(p,q) is t-independent 
and has to be evaluated in thermal equilibrium. To be explicit we use the 
self-energies ()2.2()j) in thermal equilibrium for which the matrix 5" in ()4.2j) 

reads 

S'(p,q) = J^g {cq \p~q-k)F^\k)F {cq \p) (4.3) 



1 



F (eq) G» - q - k)F {cq \k) + ^g {eq) {p -q- k)g {c<i \k) g {c(i) {p) 



with the equilibrium statistical and spectral function related by Eq. ()2.23|) . 
We have written the diagonal contribution ~ explicitly in the second 
line of Eq. ()4.2|) . since the relaxation-time approximation corresponds to 
assuming j S'SF = 0. The omission of the "off-diagonal" elements in the 
matrix S'(p, q) is a crucial feature of the relaxation-time approximation and 
responsible for its failure to describe correctly the large-time dynamics of 
thermalization as will be discussed below. 8 It corresponds to the picture that 
the dynamics of a single excitation 5F(t;p) can be described independently 
of the other modes, assuming equilibrium values for all F(q) except for the 
mode q = p. However, the true linear description is given by the full stability 
matrix S(p, q). 

Using a language of discrete momenta and a four-volume Q, for very large 
time the evolution will be dictated by the smallest eigenvalue A m i n of S/fl 
according to 

5F ~ F min e- A -* , (4.4) 

with F min the eigenvector corresponding to A m j n . The thermalization time 
can be associated with A~; n . (We omit here the complications that A min may 
be degenerate or vanish in the infinite volume limit.) The size of A m j n may be 

8 Cf. also the discussion in Ref. . 
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much smaller than the size of typical diagonal elements of the stability matrix 
or Sg Cq ' ) /2u;. The latter determines the characteristic time for a given single 
mode in the relaxation-time approximation according to Eq. (|4.7|) or (|4.10j) 
below. We also note that the eigenvector F m i n is typically not in the direction 
of the initial small deviation SF. Deviations from thermal equilibrium will 
be excited for all modes during the approach to equilibrium, explaining the 
possible failure of an approximation based on a single excitation. 

In order to compute the eigenvalues of the linearization matrix in 
Eq. ([4.2)1 . one may be tempted to diagonalize the matrix numerically from 
the regularized expression of the field-theory. 9 However, here we are only 
interested in the smallest relevant eigenvalue which determines the late-time 
behavior. In this case, it is more efficient to directly solve the time evolution 
equations (j2.15|) along the lines of Ref. ^U]. The smallest eigenvalue (or an 
upper bound for it) can then be obtained from the late-time limit of the 
dynamics. In the infinite volume limit the smallest eigenvalue may approach 
zero. The presence of vanishing eigenvalue(s) can invalidate the linear 
approach even for small perturbations around equilibrium. We emphasize 
that the direct solution of the time evolution equations (|2.15|) is not restricted 
to the linearized dynamics discussed in this section. 

Before turning to the full numerical results in Sec. EJ we consider the 
dynamics in the relaxation-time approximation. In general, the relaxation- 
time approximation does not describe the full linearized dynamics since it 
involves further assumptions. In this approximation one assumes that in the 
immediate vicinity of thermal equilibrium the dynamics is well approximated 
by omitting the contributions involving the matrix S' in Eq. (|4.2jl . We 
first discuss isotropization. Within the relaxation-time approximation the 
evolution equation for AF is given by 

relax. 

2o;-AF(t; u, q) '= E^(oj, q)AF(t; 00, q) . (4.5) 
The solution reads 

relax. 

AF{t;u,q) = c AF {u,q)e-^ iuj ' q)t , (4.6) 

9 Note that the matrix is obtained from the gradient expanded equations in the limit 
of an initial time in the remote past, i.e. to — ► —00. The latter damps out the unstable 
directions around thermal equilibrium that are present in the reversible dynamics for 
finite to. 
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where the isotropization rate is determined by the imaginary part of the 
thermal self-energy, - v '" i ' /-) •>-•-., ...,1;,,,, <.> m 



-Eg ,/2, according to 1 



clax. 

— \ time 



i( e q), 



7 iso (o;,g) * 



u,q) 



2u 



(4.7) 



We next turn to thermal equilibration. In thermal equilibrium the ratio 
of the Wigner-transformed statistical and spectral function becomes 



g^ (u,p) 



(4.8) 



according to the fluctuation-dissipation relation (|2.23jl . with the 
^-independent Bose-Einstein distribution function tibeC^)- 111 order to 
discuss the approach to thermal equilibrium, we consider the behavior of 
F j g at late times. In the relaxation-time approximation one finds 



F(t;u,p) 
g( ec (uj,p) 



relax, 
time 



n BE {uj) + - + 



Q 



Fo(u,p) 1 



-7eq(w,P)* 



• (4-9) 



Here we have used the fact that in thermal equilibrium E^/E^ = Hbe+ o- 
This yields the familiar result that the thermalization rate in the relaxation- 
time approximation is determined by the imaginary part of the thermal self- 
energy according to 



7eq(^,P) 



relax, 
time 



i(eq) 



2uo 



(4.10) 



Of course, comparing to Eq. (|4.7|) one observes the well-known result that 
the isotropization and thermalization rates agree in the relaxation-time 
approximation, i.e. 7; so = 7 eq in this case. We emphasize that in general 
isotropy is a necessary condition for thermal equilibrium while the reverse is 
not. Below we will show that beyond the relaxation-time approximation AF 
can vanish with a characteristic inverse rate, which is very different from the 
time for the approach of F to equilibrium. 



10 Note that our results fully reproduce those obtained from the linearization on the level 
of the particle number distributions. The latter can simply be obtained from (|4.2I) or (|4.5|) . 
respectively, by assuming a quasi-particle ansatz for the spectral function (i.e. a <5-function 
in Wigncr coordinates). Our equations are more general in the sense that they are not 
restricted to a quasi-particle assumption typically employed for kinetic descriptions. 
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For sufficiently weak coupling the relaxation-time rate can be computed 
from perturbation theory, which we denote by 7( pert ). For the zero- 
momentum mode this is found to be [37] 



7 



(pert) 



9g 4 T 2 



Li 2 (e- m */ T ) 



(4.11) 



where Li 2 (z) is the second poly-logarithmic function, defined by 



(z) — — / dw 
Jo 



In (1 — w) 



(4.12) 



w 



In the high-temperature limit this becomes 



7 



(pert) 



3g 4 T 2 
4:Tcm R 



(4.13) 



Our numerical results explained in the next section show that for late 
times the isotropization rate 7i SO is actually well described (on-shell) by the 
relaxation-time approximation. In particular, for sufficiently weak coupling 
the perturbative expressions for the zero-momentum mode provide a rather 
good estimate for this rate. On the other hand, thermalization occurs for 
times that exceed by several orders of magnitude the relaxation-time estimate 
()4.1L)|) . Nevertheless, examples in the literature show that for some very 
specific initial conditions the relaxation-time approximation can give a rather 
good estimate for the thermal equilibration rate as well (cf., in particular, 



One may therefore ask what are the conditions for the relaxation-time 
approximation to hold. First, one has to distinguish between quantities 
whose vanishing corresponds to an enhanced symmetry from more generic 
quantities. Once an ensemble is isotropic it always remains so at late times. 
The vanishing of AF obviously corresponds to an enhanced symmetry, i.e. 
rotation symmetry. The isotropic states may be considered as a submanifold 
of all states: Isotropization describes then the time evolution for excitations 
orthogonal to this submanifold. The relevant piece in the stability matrix S 
concerns the non-singlet representations of the rotation group. One concludes 
that isotropization is determined by the smallest eigenvalue in this submatrix 
rather than the smallest eigenvalue of 5*. (The stability matrix is block 
diagonal in the different irreducible representations of the rotation group.) 
Eigenvalues relevant for isotropization may all be much larger than the 
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smallest eigenvalue A min , which typically occurs in the singlet sector, such 
that the isotropization time is indeed much shorter than the thermalization 
time. After the system has become isotropic there still remains a slow flow 
within the isotropic subspace towards the equilibrium point. 

We next ask under what circumstances the relaxation-time approximation 
can be used to describe thermalization. Let us start with a specific initial 
condition where only one mode is excited infinitesimally. In general, due 
to the off-diagonal elements in the stability matrix S the other modes will 
be excited as well in the course of the linear evolution. This defines a 
characteristic time scale where these other modes reach amplitudes of the 
same order of magnitude as the initially perturbed mode. The relaxation- 
time approximation typically becomes invalid for larger times. On the 
other hand, for practical purposes the deviation from equilibrium may be 
already substantially reduced at that time and the following evolution of 
not much relevance. In this case the relaxation-time approximation can be 
used in practice. It has been actually observed in Refs. EES] that the 
relaxation-time approximation provides a rather good description for a single 
initially excited mode of finite amplitude, with all other modes in equilibrium. 
Since the off-diagonal elements of S' involve a volume factor Q^ 1 from the 
g-integration in Eq. (|4.2|) . the relaxation-time approximation becomes valid 
in the infinite volume limit if a set of modes of measure zero is initially 
excited. However, realistic initial conditions involve some distribution of 
initially excited modes (for example a Gaussian). Then the sum over the 
initial modes in Eq. (|4.2|) compensates, at least partially, the volume factor. 
For our initial conditions described in Sec. 12.11 even though chosen quite 
narrow in momentum space, we find that the relaxation-time approximation 
fails to describe the approach to thermal equilibrium. We believe that this 
failure happens for a very wide class of generic initial conditions. 

5 Numerical results 

In the following we present the numerical solutions of the full equations ()2.15|) 
with the self-energies (|2.20|) . The results are then confronted with the 
estimates from the relaxation-time approximation discussed above. We 
employ a rather weak coupling g = 0.5 and consider a wide range of 
anisotropic initial conditions (|2.12j) - (J2.13j) with amplitudes A = 10 — 160. 
The width a of the initial particle number distribution no(p) of Eq. (|2.13|) 
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Figure 5: Isotropization. On the left the normalized anisotropy correlator 
AF(t,t;pt s ) is shown as a function of time for different initial amplitudes A with 
same energy density. One observes an exponential decay of the anisotropy with 
AF — > at late times. While the behavior at early times is sensitive to the initial 
amplitude, the isotropization rates become universal at late times. On the right 
the corresponding pressure differences normalized to the energy density are shown. 



is chosen such that the energy expectation value is the same for the 
different initial conditions. Therefore, if the system thermalizes then the 
late-time results are determined by the same temperature. The employed 
energy density corresponds to that of a thermal system with temperature 
T = 1.03 mij, where we always express dimensionful quantities in units of 
the renormalized thermal mass m#. 

We first consider the anisotropy correlator AF(t±, t2!Pts) defined in 
Eq. (12.16}) with p ts = 0.93 rriR. At equal times t\ = £2 it corresponds to 
the frequency integral of the corresponding correlator in Wigner space, 
i.e. AF(t,t;p ts ) = J[du/2n} AF(t; 0J,p ts ) (cf. Sec. EJ). The left graph of 
Fig. El shows the time evolution of AF(t,t;p ts ) normalized to 2F(t,t;p ts ) 
for different anisotropy amplitudes A. One observes that at early times the 
behavior of AF/2F depends on the initial amplitude. However, for times 
t > 300/m^ the isotropization rate becomes independent of the details 
of the initial conditions. It approaches an exponential behavior with rate 
7iso(pts) = 0.0061(2) m R . 

A similar behavior can be seen in the right graph of Fig. E] for the 
pressure anisotropy AP(t) normalized to the energy density e. Its late- 
time isotropization rate turns out to be rather accurately described by the 
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characteristic rate of the maximally populated mode, i.e. 7is (pts)- The 
pressure anisotropy and energy density have been calculated form the 
diagonal components of the energy-momentum tensor (j2.21|) as e = T 00 
and AP = T 33 — T 11 = T 33 — T 22 . A value reasonable for the use of 
hydrodynamics, say AP/e < 0.1, is reached particularly fast for large A. 

In order to directly compare the results to the relaxation-time 
approximation (J4.7j) . we consider the late-time evolution of the Wigner- 
space correlator AF (t; u , p ts ) ■ Following the results of Sec. EJ AF(t;uj,p ts ) 
is described by the evolution equation (|3.3|) for sufficiently late times. If 
the relaxation-time approximation is valid then Eq. ()3.3|) can be replaced 
by the simpler expression of Eq. (J4.5|) . Comparing the two equations one 
observes that in this case the ratio — (EgAF+AEgF— E^A^— AY> F ~p)/2u:AF 
must approach the relaxation-time rate £g CC ^/2u/ eq ). To establish the 
equivalence, we first evaluate the former expression using the full solution 
of the nonequilibrium evolution from Eqs. (j2.15j) and (|2.2()jl . In a second 
calculation, we evaluate Eg /2u/ eq ) directly in thermal equilibrium for the 
same energy density. The results of these calculations are shown in Fig. Q 
for on-shell u, i.e. evaluated from the peak of the spectral function JT] 
for momentum p ts . The lower curve of Fig. shows the nonequilibrium 
time evolution of the ratio — (S^AF + AT, F — S^A^ — AHpQ)/2uAF. 
One observes that indeed for times t > 500/m^ the nonequilibrium ratio 
approaches closely — T^f"^ /2u ( - eq \ where the latter value is indicated by the 
square symbol on the right of the figure. This verifies standard relaxation- 
time approximations for the description of isotropization at sufficiently late 
times. 

It should be stressed that this agreement of the nonequilibrium 
isotropization rate with the normalized imaginary part of the equilibrium 
self energy, i.e. — /2u^ eq \ is very nontrivial for relatively early times 
t > 500/m/j: The imaginary part of the nonequilibrium self-energy at that 
time is still far from equilibrium! This is demonstrated in Fig. ^ where 
the upper curve shows the normalized imaginary part of the nonequilibrium 
self-energy, or —H e (t',u,q)/2u(t). The calculation (solid line) was performed 
until t = 50000/rnR, and the dotted line in the figure represents a fit. 11 
Only at very late times, i.e. much later than the inverse relaxation-time rate 

n The plotted fit is obtained using a power-law behavior, however, the available late- 
time range of numerical data can be fitted using an exponential behavior with comparable 
accuracy. 
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Figure 6: Equal-time correlator F(t,t;p) = F(t,t;pj_,P\\) for different three- 
momenta as a function of time in units of the renormalized thermal mass rriR. The 
lines show the modes for transverse, F(t,t;p±,0), and for longitudinal momenta, 
F(t, t; 0,pi|), for the values of p± and pn as indicated in the figure. One observes 
isotropy with F(t,t;p±,0) = F(t,t;0,p\\) for p± = pn to good accuracy for 
times t > 300/rriR. In contrast, equilibrium has not been achieved at that 
time and the effectively isotropic modes continue to evolve towards a thermal 
distribution. For comparison, the dashed line of the inset shows the Bose-Einstein 
(BE) distribution along with the nonthermal time-dependent distribution function 
n(\k\) as a function of \k\/m,R at time t = 300/mR. 

t ^> — 2co> < - cq ' l /X!g eq \ its nonequilibrium evolution comes rather close to its 
thermal value. This points out that the relaxation-time approximation fails 
to give an estimate for the thermalization time. According to the relaxation- 
time approximation discussed in Sec.EJ the characteristic thermalization time 
agrees with the isotropization time — 2u/ eq ) /T,^ q \ This is sharp contrast with 
the findings represented in Fig.^showing a large separation of times between 
isotropization and thermalization. We emphasize that for the considered class 
of initial conditions the large separation is also found if the coupling is further 
weakened (cf. also the discussion in Sec. HJ). 

The faster isotropization as compared to thermalization can also be 
observed from the real-time correlator modes F(t,t;p) in Fig. |3 Isotropy 
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Figure 7: The left graph shows the evolution of the amplitudes and the right 
graph shows the damping rates extracted from the unequal-time correlators F and 
p as a function of t, as defined in Eq. 1)5.1 Jl for the maximally amplified mode with 
\p\ = pts and for a low-momentum mode. 



implies F(t, t;p±,0) = F(t,t;0,pu) for p± = py, which is achieved to good 
accuracy for times t > 300/m#. In contrast, the modes still show strong 
deviations from thermal equilibrium at that time. To emphasize this fact 
the inset compares the time-dependent nonthermal effective particle number 
distribution [SSI EEU Clj n(k) — n(\k\) as a function of \k\/m,R at time 
t = 300/m^ with the Bose-Einstein distribution. 

The amplitudes and damping rates plotted in Fig.[7|are the fit parameters 
Cf/p and jF/ p (t;p), respectively, in 

F{t u t 2 ;P) = C? f (t;p>Mi;p)(fi-t 2 )]e^^( f1 ^ 
P {h,t 2 )V) = Cp^^sin^j^^-^le-^^-* 3 ) (5.1) 

The values are extracted at each t — (t\ + t 2 )/2 adjusting Cp/ p {t;p), 
Up/ p (t',p) and ■JF/ P (t;p). We observe (not shown in the plot) an isotropic 
Up/p(t',p) parameter form the very beginning. In thermal equilibrium the 
damping parameters of Eq. (|5.1|) would equal half the rate — £ e (i; lu)/2lu if 
the correlator modes F(t; lu, p) and g(t; lu, p) were of exact Breit-Wigner form. 

6 Conclusions 

Isotropization happens much faster than the approach to thermal 
equilibrium. This is not a surprising fact, since only a subclass of those 
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processes that lead to thermalization are actually required for isotropization: 
"ordinary" 2 «-> 2 scattering processes are sufficient in order to isotropize 
a system, while global particle number changing processes are crucial to 
approach thermal equilibrium. For the massive weak-coupling scalar theory, 
the 2^2 processes are of order g 4 , while global number changing processes 
such as f <-» 3 processes are of order g 8 . For small couplings the 2^2 
scatterings can, therefore, dominate the early-time behavior. While these 
processes are in principle sufficient to achieve complete isotropization, they 
are not sufficient to provide a quantitative description of thermalization. 
In the absence of global particle number changing processes they would, 
for instance, lead to a spurious chemical potential in the Bose-Einstein 
distribution at late times, which is clearly absent for the thermal equilibrium 
theory of real scalar fields. 

We demonstrate that these apparently small higher order contributions 
lead to quantitatively important corrections for general nonequilibrium 
situations. To establish this requires to go beyond the relaxation-time 
approximation, which fails to describe the dynamics close to equilibrium 
in general. We emphasize in Sec. 0] that the relaxation-time ansatz neglects, 
in particular, all off-diagonal elements in the stability matrix around thermal 
equilibrium. Whether the true small eigenvalues of the stability matrix play 
an important role for the late-time behavior depends on the physical quantity 
of interest. We have shown that the isotropization rate is at sufficiently 
late times well described by the relaxation-time approximation, though the 
latter fails even to give a correct order-of-magnitude estimate for the thermal 
equilibration time. 

Since our approach does not involve a gradient expansion, we can 
determine the range of validity of a description to lowest-order in gradients, 
which is typically employed in kinetic equations and to motivate the 
relaxation-time ansatz close to equilibrium. The earliest time for the 
applicability of the lowest-order gradient expansion is found to be of the 
order of the characteristic isotropization time t iso . Therefore, the gradient 
expansion and, consequently, the respective Boltzmann equation cannot be 
used to compute the isotropization time. Similar problems with a lowest- 
order gradient expansion concern thermalization. 

Far-from-equilibrium isotropization can be very important for the 
understanding of collision experiments of heavy nuclei. An isotropic equation 
of state P/e is a crucial ingredient for the apparently successful application 
of hydrodynamic descriptions. We have shown in a previous publication 
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that the prethermalization pQ of the equation of state to its equilibrium 
value is insensitive to the thermal equilibration time r eq . The relevant 
time scale for the early validity of hydrodynamics could then be set by the 
isotropization time, which can be much shorter than r eq . In this context 
it is interesting to note that our simple scalar theory indicates that large 
initial anisotropics can lead to enhanced characteristic isotropization rates 
as compared to conventional relaxation-time rates. 
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